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Abstract 

The contributions of complex rays and the secondary radiation 
shed by surface waves to scattering by a dielectric sphere are 
calculated in the context of the Debye series expansion of the Mie 
scattering amplitudes. Also, the contributions of geometrical rays 
are reviewed and compared with the Debye series. Interference 
effects between surface waves, complex rays, and geometrical rays are 
calculated and the possibility of observing these interference 
effects is discussed. Experimental data supporting the observation 
of a surface wave - geometrical ray interference pattern is 
presented. 


Introduction 

The infinite series of transverse electric (TE) and transverse 
magnetic (TM) spherical multipole partial waves, known as the Mie 
scattering formalism, is an exact solution to the scattering of a 
linearly polarized plane electromagnetic wave by a dielectric 
sphere . 1 ' 3 Being an exact solution, the Mie infinite series contains 
all the effects that contribute to the scattering. These effects for 
the most part, are not readily identifiable in the complicated terms 
of the Mie infinite series. It turns out that writing each term of 
the Mie infinite series as another infinite series, known as the 
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Debye series, clarifies the physical origins of many effects that 
occur in electromagnetic scattering. 4 ' 9 In ray theory, when a 
geometrical light ray is incident on a dielectric sphere, it is 
partially reflected by the sphere surface, partially transmitted 
through the sphere, and partially transmitted after making an 
arbitrary number of internal reflections. Analogously, each term of 
the Debye series decomposition of an individual TE or TM partial wave 
scattering amplitude may be interpreted as diffraction of the 
corresponding spherical multipole wave, or its reflection by the 
sphere surface (p = 0) , or transmission through the sphere (p = 1) , 
or transmission after making p-1 internal reflections (p > 2) . 

Expressed in this way, the total scattered electric field takes 
the form of a double sum. One sum is over partial waves, and the 
other sum is over the number of interactions p that each partial wave 
makes with the sphere surface before propagating into the far-field. 
For scattering in the short wavelength limit, the sum over partial 
waves has long been known to be slowly convergent. But the sum over 
surface interactions is expected on physical grounds to be rapidly 
convergent for nearly all partial waves since summing over multiple 
internal reflections resembles the summing of a geometric series. 
The major exception to this rapid convergence is for partial waves in 
the edge region, corresponding to geometrical light rays incident on 
the sphere at grazing incidence, where the internal reflection 
coefficient is nearly unity. 5 

For most scattering angles, the far-field scattered electric 
field is dominated by the contributions corresponding to geometrical 
light rays. 10 Occasionally, however, other mechanisms dominate the 
scattering. For example, the p-1 order rainbow caustic occurs at the 
scattering angle where two light rays coalesce after having made p-1 
internal reflections within the dielectric sphere. 11 A rainbow 
separates two regions of scattering in which the number of 
contributing geometrical light rays differs by two. To one side of 
the rainbow, the two p-1 internal reflection rays form an 
interference pattern known as supernumerary rainbows. 12 Although 
these two rays are absent on the other side of the rainbow, their 
contribution to the scattered light does not discontinuously fall to 
zero at the two-ray to zero-ray transition. Rather, their 
contribution smoothly but rapidly decreases due to their 
metamorphosis into a "complex ray" in the zero-geometrical ray 
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region . 

Another example where the scattering is dominated by effects 
other than geometrical light rays occurs when light is incident on 
the sphere at grazing incidence. Beyond the scattering angle of the 
grazing incidence ray, the contribution of this ray to the scattered 
intensity again does not discontinuously fall to zero at the one-ray 
to zero-ray transition. This is due to the creation of 
electromagnetic surface waves at the point of grazing incidence on 
the sphere . 5 ' 7,16 ' 17 These waves travel along the sphere circumference, 
and while doing so, they shed secondary radiation which propagates 
into the far-zone. This secondary radiation produces radiation 
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damping of the surface waves. As a result, the surface wave 
amplitude dies off exponentially along the sphere's circumference and 
the amplitude of the secondary radiation propagating into the far- 
zone dies off exponentially as a function of scattering angle. The 
far-field intensity produced by the surface waves joins smoothly with 
the far-field intensity produced by the geometrical light rays in the 
vicinity of the one-ray to zero-ray transition. This smooth 
intensity transition is known as the Fock transition. 5 " 7,18 

A question of interest is whether under favorable circumstances, 
some of these non-geometrical mechanisms that weakly contribute to 
the scattering, namely surface waves and complex rays, may be 
observed. It has been known for some time that the optical glory is 
dominated by the contribution of p = 2 surface waves. 7 ' 15,19 ' 20 In 
addition, a subtle effect of the complex ray in the neighborhood of 
a transverse cusp caustic has been predicted 21 " 24 and tentatively 
observed. 25 One of the purposes of this paper is to assess the 
practicality of additional observations of surface waves and complex 
rays in light scattering by a dielectric sphere. The observations 
that we examine employ the interference of surface waves or complex 
rays with geometrical light rays. Observing relatively weak 
contributions to scattering via their interference with more dominant 
contributions has long been employed in quantum mechanical 
scattering. 26 We believe that this method has not been previously 
considered for surface waves or complex rays in the context of light 
scattering by a dielectric sphere. 

This paper has a second and more theoretical purpose as well. 
Although light scattering experiments measure the full scattered 
intensity rather than only one Debye series component at a time, the 
Debye series decomposition of the scattering amplitudes is a powerful 
tool for understanding the physical mechanisms that produce the 
scattering. The individual Debye component intensities allow one to 
examine a single scattering mechanism in isolation from all the other 
mechanisms that either dominate it or otherwise obscure its effects. 
The Debye series component intensities are found to exhibit a number 
of novel features and interference structures that are not evident in 
the total Mie intensity. Even though many of these novel structures 
may not be observable in practice, we believe that determining their 
physical origins increases our fundamental understanding of the 
scattering process. 

The remainder of this paper is organized as follows. In 
section 2 we review the contribution of geometrical light rays to the 
scattered intensity. In section 3 we review the Airy theory of the 
p-1 order rainbow and compare its accuracy to the p-term of the Debye 
series expansion of the scattered intensity. We parameterize the 
complex ray in the zero-ray region adjacent to the rainbow in terms 
of the Airy integral. In section 4, we review the complex angular 
momentum parameterization of surface waves and compare its accuracy 
to the p = 1 (i.e. transmission) term of the Debye series. Next, in 
sections 5-7 using these models for complex rays, surface waves, and 
geometrical light rays, we examine the interference between 1) 
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surface waves and geometrical rays; 2) complex rays and geometrical 
rays; and 3) surface waves and complex rays that occur in the various 
Debye series component intensities. In section 8, we describe an 
experiment in which we observed the surface wave - geometrical ray 
interference. Finally, we present our conclusions in section 9. 

Geometrical Light Rays 


Consider a linearly polarized geometrical light ray with field 
strength E 0 and wavelength X incident with the angle 6 - t on a 
dielectric sphere of radius a and refractive index n as in Fig. 1. 
The size parameter of the sphere is defined as 


x = 


2 ma 


( 1 ) 


The angle of deflection of the geometrical ray as it leaves the 
sphere after having made p-1 internal reflections is 


0 = (p-l)m + 20 i - 2p0 r 
where 

sin0 i = n sin0 r . 

The scattering angle corresponding to this deflection is 


0 


■( 


0 - 2m N if 2nN <. 0 £ 2m (N+V 2 ) \ 

2m (N+l) - 0 if 2 k(N+Vz) <0£2m(W+l)j 


( 2 ) 


(3) 


(4) 


where N is an integer. This relation confines the scattering angle 
to the interval 0 < 0 < 180°. The fraction of the geometrical ray's 
incident TE or TM polarized electric field that is transmitted from 
the exterior to the interior of the sphere is T 21 (00, the fraction 
internally reflected is R n (0j), and the fraction transmitted from the 
interior to the exterior is T I2 (0j). The superscript 1 denotes the 
region inside the sphere and the superscript 2 denotes the region 
outside the sphere. When x » 1, these fractions are approximated by 
the Fresnel coefficients for obligue incidence on a flat interface 
with either the TE or TM polarization. 27 

The optical path length of the ray from the entrance plane of 
the sphere to the exit plane is 


L = 2a(p n cos0 r - cos0 i ) + 2a. 


(5) 


The electric field of the ray scattered in the 6 direction after p-1 
internal reflections is then" 


E n a 


T?P - 0 

geometrical ray \ v ' 


sin0 i cos0 i 

2sin0 

p COS0,- 

1 - 

n cos0 r 


t 21 (0^) [i? 11 ^.)]*- 1 r 12 (0^.) 


x gikR ^2niL/\. gi C 


( 6 ) 

for either polarization state where R is the distance from the center 
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FIGURE 1: The deflection of a geometrical light ray through the angle 
9 by a dielectric sphere of radius a and refractive 
index n. 
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of the sphere to the observer and f is a phase factor that contains 
the effect of the ray trajectory crossing focal lines. 28 For almost 
all values of p, there are a number of values of 9 , that give rise to 
the same scattering angle 9 . As a result, when eq. (6) is used in 
obtaining the total ray optics scattered electric field, all of the 
contributing ray trajectories for the scattering angle 6 must be 
summed . 

As mentioned earlier, the far-field Mie scattered electric field 
at most angles is dominated by the contributions of geometrical light 
rays. A comparison between the geometrical ray model and the various 
Debye series component intensities can be used to map the angular 
regions corresponding to the various p-terms that are dominated by 
the geometrical light rays and the angular regions that are dominated 
by other effects. The appropriate formulas for the various terms of 
the Debye series expansion of the scattered electric field and their 
relation to the TE and TM polarization states of the geometrical rays 
are given in Appendix A. 

Figure 2a-d shows the comparison between the Debye scattered 
intensities (Il=|Sl| 2 ) and the corresponding TE polarization state of 
the geometrical ray for p = 0, 1, 2, and 3 for a sphere with x = 100 
and n = 1.333. Note in Fig. 2a that the Debye series plot for p = 0 
also includes the diffraction term. Details describing the necessity 
of combining the p = 0 [reflection] term with the diffraction term 
are given in Appendix A. Also note in Fig. 2a that agreement between 
the geometrical ray model of reflection and the Debye p = 0 plus 
diffraction intensity is excellent for the larger scattering angles 
where the reflection term of the Debye series dominates. Thus, ray 
optics is a good model for reflection from a sphere with size 
parameters as low as 100 or less. 29 

Figure 2b compares the geometrical ray model of the transmitted 
light with the p = 1 term of the Debye series. The ray optics model 
predicts that no light is scattered past 9 = 82.79° whereas the 
intensity contributed by the p = 1 term of the Debye series continues 
out to 180°. This continuation is produced by the secondary 
radiation shed by surface waves and will be addressed in more detail 
in section 4. One other point of interest is the oscillations near 
9 = 180° in the Debye series component intensity. This is the 
interference pattern that surrounds the backscatter glory axial 
caustic. 30 " 32 

Figures 2c-d compare the geometrical ray model of the light 
making one and two internal reflections with the p = 2 and the p = 3 
terms of the Debye series. In both cases ray optics predicts large 
angular intervals where no rays are scattered and predicts an 
infinite intensity at the positions of the rainbows. In contrast, 
the Debye series component intensities extend over all scattering 
angles and show a broad increase in intensity near the rainbow with 
the characteristic supernumeraries located to one side. Also note in 
Fig. 2c the oscillations in the intensity near 9 = 0 ° for the p = 2 
Debye contribution. This is the interference pattern that surrounds 
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a) Debye p = 0 plus diffraction 
compared with the reflected 
geometrical rays; 



c) Debye p = 2 compared with the 
single internal reflection 
geometrical rays; 



b) Debye p = l compared with the 
transmitted geometrical rays; 



d) Debye p = 3 compared with the 
two internal reflection 
geometrical rays. 


FIGURE 2: Comparison of the p-term Debye series component intensity 
Ij with the prediction of geometrical optics, eq. (6) for 
all the contributing TE polarized ray trajectories for 
x = 100 and n = 1.333. 
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the forward glory axial caustic. 26,33 


Airy Theory of the Rainbow 


As is seen in Figs. 2c-d, the electric field of eq. (6) for the 
geometrical rays that have made p-1 internal reflections incorrectly 
predicts an infinite intensity in the direction of the rainbow 
scattering angle 6 R . This angle is given by eq. (4) with 

0* = (p-1) n + 2 0? - 2 p 0* (7) 


and 


cos 2 0? = 


n 2 

P 2 


1 

1 ' 


( 8 ) 


In Airy theory, the shape of the wavefront leaving the dielectric 
sphere in the vicinity of 6 R is cubic to a first approximation. 34 If 
this exiting wavefront is then Fraunhofer diffracted into the far- 
zone, the TE or TM polarized electric field in the vicinity of 9 R 
becomes 35 


-'Airy 


( 0 ) 


E 0 a e 


ikR 


2 ix sin0? 
sin0* 


*1 r 21 (0?) 

h 3 


[j? 9 * 11 (Of) ] p_1 r 12 (0?) 


where 

A = 0 - 0 R , 


x Ai 


2 \ 
-x 3 A 


i 



e 2TtiL R /X 



(9) 


( 10 ) 


h = 


(p 2 -!) 2 

P 2 


(p 2 ~P 2 ) 1/2 
(n 2 -l ) 3/2 ' 


( 11 ) 


L r is eq. (5) evaluated at the rainbow angle, and Ai is the Airy 
integral. 36 For negative values of its argument, the Airy integral is 
oscillatory and describes the supernumerary interference pattern in 
the two-geometrical ray region. For positive values of its argument, 
the Airy integral rapidly dies off in the zero-geometrical ray region 
as 


Ai{u) ~ 



for u >> 1. 


( 12 ) 


In the complex angular momentum analysis of the Debye series 
terms, for a given value of p, the sum over partial waves is 
converted into an integral over an effective impact parameter 
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(i.e. how far off center an incident ray strikes the sphere) via the 
modified Watson transformation. 5 ' 7 The numerical value of this 
integral is dominated in the small wavelength limit by the regions of 
stationary phase and residue poles of the integrand. The stationary 
phase regions produce contributions to the integral that resemble the 
effects of geometrical light rays. This result is the motivation for 
the localization principle mentioned in the next section. The 
residue pole contributions correspond to the shedding of secondary 
radiation by surface waves. 

After the two geometrical rays coalesce at the rainbow 
scattering angle, if the scattering angle is further decreased, the 
impact parameters of the two rays in the complex angular momentum 
analysis leave the real axis and migrate into the complex plane. The 
contribution to the Debye term integral from these complex impact 
parameter stationary phase points is known as the complex ray 
contribution to the scattering. 

The complex ray is loosely analogous to an exponentially damped 
wave or evanescent wave. When a plane wave is incident on a flat 
dielectric interface, a portion of the wave is refracted through the 
interface at an angle 9 r as determined by Snell's law 
(0 r = arcsin (n ; /n r sin 9,). If the wave is incident from a denser 
medium (i.e. nj/n r > 1) and angle of incidence 9, of the wave increases 
past the critical angle for total internal reflection 
(i.e. n ; /n r sin 9, > 1), then 9 r must become complex to satisfy Snell's 
law. Thus the angle of refraction leaves the real axis and migrates 
into the complex plane and an evanescent wave is formed. The 
imaginary part of the angle of refraction is responsible for the 
damping of the evanescent wave. 

Since the zeroth order approximation to the complex angular 
momentum analysis gives results identical to Airy theory, we take the 
rapid falloff of the Airy integral for u » 1 as our parameterization 
of the scattering produced by the complex ray. 

The intensity of the p-1 order rainbow in Airy theory 
ifiryO) =| E% ily (Q)\ 2 ( 13 ) 

has been previously compared with the full Mie intensity for 
x - 15,000 for only the first and second order rainbows (p = 2, 3). 35 
The reason that comparisons for higher order rainbows have not been 
previously made is that the Mie intensity is dominated by the much 
larger diffracted, reflected, and transmitted contributions in the 
regions where these rainbows occur. Only in Alexander's dark band 
between the first and second order rainbows is the background 
intensity low enough so that the features of the fifth and sixth 
order rainbows may be qualitatively examined using Airy theory. 35,37 ' 38 
However, if I p Airy (0) were compared with the intensity corresponding to 
only the p-term of the Debye series expansion of the Mie scattering 
amplitudes, the other contributions to the Mie amplitudes that would 
normally obscure the relatively weak high order rainbows would be 
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removed from consideration. In the remainder of this section we use 
the Debye term comparison to provide a sensitive test of the validity 
of Airy theory for high order rainbows. 

Using eqs.(9,13) and the results of Appendix A, the comparisons 
between Airy theory and the Debye series terms are given for the 
first through fifth order rainbows in Figs. 3-7 for x = 1000 and 
n = 1.333 for the dominant polarization (i.e. for the Debye 
scattering intensity 1,(0) and the TE polarization for Airy theory). 

The Airy theory formula of eq. (9) neglects the variations in the 
Fresnel reflection and transmission coefficients as function of 9, in 
the vicinity of 9*. This variation is relatively unimportant for the 
primary rainbow in the dominant polarization as is seen in Fig. 3. 
But, as evidenced by Figs. 4-7, the variation in the Fresnel 
coefficients is important for all higher order rainbows since 0 R , is 
closer to n/2 where the transmission and reflection Fresnel 
coefficients are rapidly increasing or decreasing functions of 0 ,. 
This variation is also important in the nondominant polarization 
(i.e. for the Debye scattering intensity I 2 (0) and the TM polarization 
for Airy theory) since 0 R is near the Brewser angle for this 
polarization and the internal reflection of the contributing rays is 
weak. 39 For very large droplets with x « 20,000 for the p-1 order 
rainbow in the dominant polarization, the Airy theory approximation 
of eq.(13) was found to be virtually identical near the primary 
maximum to the exact Debye series calculation for all the values of 
p that were examined. 


Surface Waves 

For x >> 1, the localization principle 40 associates a small 
number of spherical multipole partial waves centered about the 
partial wave number £, vc with a geometrical light ray whose angle of 
incidence on the sphere is 

sine. * 1™. (14) 

1 x 

For example, small values of £ correspond to rays incident near the 
centerline of the sphere and values of £ - x correspond to light rays 
incident on the sphere at grazing incidence. In ray theory, such a 
grazing ray is entirely reflected by the sphere. In wave theory, the 
grazing incidence also creates surface waves which propagate along 
the circumference of the sphere. At every point along the 
circumference, they shed secondary radiation tangentially which then 
propagates into the far-zone. A portion of the surface waves also 
refract into the sphere at the critical angle. The angle of the 
refracted wave inside the sphere is given by Snell's law. This angle 
turns out to be the critical angle because the incident wave (i.e. 
the surface wave) is propagating parallel to the surface. After the 
refracted wave propagates across the interior of the sphere and 
reaches the opposite side, it is again refracted at the surface. 
Since the angle the refracted wave makes with the surface is again 
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FIGURE 3: Comparison of the p = 2 Debye component intensity I, with 
the Airy theory approximation of eq.(13) for the TE 
polarization state, x = 1000, and n = 1.333. 
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FIGURE 4: Comparison of the p = 3 Debye component intensity I, with 
the Airy theory approximation of eq. (13) for the TE 
polarization state, x = 1000, and n = 1.333. 
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FIGURE 5: Comparison of the p = 4 Debye component intensity I, with 
the Airy theory approximation of eq. (13) for the TE 
polarization state, x = 1000, and n = 1.333. 
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FIGURE 6: Comparison of the p = 5 Debye component intensity I, with 
the Airy theory approximation of eq. (13) for the TE 
polarization state, x = 1000, and n = 1.333. 
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FIGURE 7: Comparison of the p = 6 Debye component intensity I, with 
the Airy theory approximation of eq.(13) for the TE 
polarization state, x = 1000, and n = 1.333. 
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the critical angle, it exits the sphere parallel to the surface and 
forms another surface wave. Thus a portion of the original surface 
wave takes a shortcut through the sphere before returning to the 
surface again to shed more secondary radiation tangentially into the 
f ar-zone . 3,6,18 

The phenomenon of surface waves has long been known and has many 
applications in electromagnetism, 13,41 acoustics, 42 and quantum 
mechanical scattering. 43 It is only comparatively recently, 7 however, 
that the scattering amplitude for electromagnetic surface waves on a 
dielectric sphere has been calculated. In this section we review 
these results and compare them with Debye series calculations. 


The surface waves, otherwise known as creeping waves or surface 
guided modes, resemble electromagnetic fields propagating in a leaky 
waveguide duct on the exterior of the sphere. 44 The energy which 
leaks out of the effective waveguide duct into the sphere is the 
previously mentioned "shortcut through the sphere". The energy which 
leaks out of the sphere surface and propagates into the far-zone is 
the "shed secondary radiation" in this section. We consider the 
surface wave in the j waveguide mode, where j > 1 is an integer, that 
takes p shortcuts through the sphere and thereby makes p-1 internal 
reflections. Such a surface wave is pictorially represented in 
Fig. 8. Let £ be the total angle traveled by the surface wave on the 
circumference of the sphere. The distance 

T = — (n 2 -l) 1/2 (15) 

n 

is the length of a shortcut through the sphere made at the critical 
refraction angle 6 r c given by 


sin0£ = — • 
n 


^ shortcut 


= nk, 


L rJ 

^ surface 


= k + (*\ 1/3 + i ■ 3 - 1 - (-) 1/3 - - (n 2 -l) 1/2 

2 a \ 2 / [ 2a \2) a J 


where Xj is defined by 
Ai (-Xj) = 0 
and 





(16) 

the 

shortcut is 





(17) 

the 

circumference 

is 7 

(n 2 - 

1)1/2 


(18) 




(19) 


1 for the 5. scattering amplitude ) 

K = I . (2°) 

n 2 for the S 2 scattering amplitude J 

The real part of k j surf , ce indicates that the propagation speed of the 
surface wave along the circumference is less than c. The imaginary 
part describes the radiation damping. The scattered electric field 
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FIGURE 8: Secondary radiation shed into the far-zone by a p = 2 
surface wave that travels the angular distance £ = £, + £ 2 
along the circumference of the sphere. The path segments 
marked T denote the shortcuts made by the surface wave 
through the sphere and the thick arc segments denote 
propagation along the sphere surface. 
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due to the secondary radiation produced by the j mode of the p-1 
internal reflection surface wave is 7 


rpP. J 

^ surface wave 


( 0 ) 


E 0 a e 


ikR 




i?(sin0) 1/2 2n 1/2 aj 


?(i) 


2 V 1 / 6 


Q^-P^shorccuc^ Q^^suzface a ^ q 2 1 k a 


( 21 ) 


(p-1) ! 


m= 1 


(m-1) ! ( p-m ) 


2k 


(n 2 -i) V2 


ml 


where a'j is defined in terms of the derivative of the Airy integral 
by 


Ai'i-Xj) = a/ (22) 

and the relationship between the scattering angle and the angle £ 
along the circumference is given by eq. (4) with 

0 = pn -2 p0£ + £. (23) 

The associated scattered intensity for either the S, or S 2 scattering 
amplitude is 


I 


p 

surface wave 


( 0 ) 


E surface wave ( 0 ) 


7*1 


(24) 


In practice, only the j = 1 term of this sum is usually considered. 
This corresponds to the largest amplitude surface wave and the 
strongest secondary radiation in the far-zone. 


The accuracy of the surface wave parameterization of eq. (21) was 
tested by comparing 1) the scattered intensities due to both the TE 
and TM polarized p = 1 geometrical rays plus their corresponding 
p = 1, j = 1 surface waves to; 2) the intensity of the p = 1 term of 
the Debye series for the Sj and S 2 scattering amplitudes respectively. 
Figures 9a-b show the results for both polarization states for 
x = 1000 and n = 1.333. In the vicinity of the critical scattering 
angle, a complete comparison would require that eqs.(6,21) be 
replaced by the Fock transition formulas of Refs. 5-7 which smoothly 
connect the ray theory intensity below the critical scattering angle 
with the surface wave intensity above the critical scattering angle. 
As x increases, the surface wave encroachment of the scattered 
intensity into the scattering angle region 9 > 9 C decreases. 


Surface Wave - Geometrical Ray Interface 

Consider the vicinity of the p-1 order rainbow. In ray theory, 
the geometrical light ray that is incident on the sphere with the 
angle 9 R exits at the rainbow scattering angle 0 R . The two rays that 
exit in the direction 9 > 9 R and interfere to produce the 
supernumerary pattern are incident on the sphere with the angles 
9 j < 9 R | and 9, > 9 R , respectively. As one progresses further into the 
supernumerary region, the angles of incidence of the contributing 
rays continue to decrease and increase respectively until the one 
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FIGURE 9: Comparison of the p = 1 Debye component intensity with the 
transmitted geometrical ray and surface wave contributions 
of eqs.(6,21) for x = 1000 and n = 1.333. 

a) the I, Debye scattered intensity, the TE polarized ray, 
and the TE polarized surface wave; 

b) the I 2 Debye scattered intensity, the TM polarized ray, 
and the TM polarized surface wave. 
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with the larger incident angle approaches 9 { = n/2. In ray theory, 
the supernumerary interference pattern ends at 9, = 77/2. This 
corresponds to the scattering angle 9 C given by eq.(4) with 

0 C = pn - 2 pQ c r . (25) 

For scattering angles beyond 9 C only a single remaining geometrical 
ray, the one with 0-, < 9 R ,, contributes to the scattered intensity. 
In wave theory however, the grazing incidence geometrical ray 
launches surface waves. For 9 > 9 C , the remaining geometrical ray 
and the secondary radiation shed by the surface waves continue to 
interfere, thereby extending the supernumerary region beyond the ray 
theory limit. This surface wave - geometrical ray interference is 
indicated pictorially in Fig. 10. 


For the first order rainbow with n = 1.333, the supernumerary 
region begins at 9 R = 137.92°, and ends at 9 C = 165.57° according to 
ray theory. The interference of the p = 2 surface waves with the 
remaining p = 2 geometrical ray having 9, < 14.6° encroaches into the 
glory region and is not easily observed because of its interference 
with other backscatter ing mechanisms. For the third order rainbow, 
the supernumerary region begins at 9 R = 41.74° and ends at 9 C = 28.86° 
according to ray theory. The interference between the p = 4 surface 
waves and the remaining p = 4 geometrical ray having 9, < 62.6° 
encroaches into the forward diffraction region and is again not 
easily observed because of the dominance of other forward scattering 
mechanisms. For the second order rainbow however, the effect may be 
amendable to observation. For p = 3 the supernumerary region begins 
at 9 R = 129.11° and ends at 9 C = 111.64° according to ray theory. The 
surface wave - geometrical ray interference propagates toward 
9 = 90°, where the total Mie intensity, and thus the strength of the 
competing scattering mechanisms is generally low, and the effect 
might be observable. 

In order to see how the interference is imprinted on the 
scattered intensity, we compare in Figs, lla-b, the p = 3 term of the 
Debye series for each polarization with 


I approximate (0) I -E geometrical ray ( 0 ) 


„p=3, j=l 

s -‘ surface wave ' w > 


(26) 


using only one of the geometrical rays that contribute to the second 
order rainbow, the one incident closer to the centerline of the 
sphere. In this figure, eq. (26) has been offset by a factor of 100 
for clarity. As mentioned earlier, the contribution of the 
geometrical ray becomes infinite at the rainbow angle. This is 
apparent in Figs. 11a, b. The Debye series component intensity in 
Fig. 11 shows an oscillatory structure for both polarizations that 
extends well below 9 C . This is the surface wave-geometric ray 
interference. The approximation of eq. (26) agrees well with the 
Debye series result for 9 « 9 C both in the amplitude and period of 
the oscillation. Near 9 C the comparison is poor because eq. (26) 
should be replaced by the previously mentioned Fock transition 
formulas at the geometrical ray-to-surf ace wave transition. 5 The 
falloff in the amplitude of the oscillatory structure of the Debye 
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FIGURE 10: A p = 3 geometrical ray (solid line) and a p = 3 surface 
wave (dashed line) that interfere in the far-field. 
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FIGURE 11: Comparison of the p = 3 Debye component intensity with the 
approximation of eq. (26) for x = 1000 and n = 1.333. The 
approximation of eq. (26) has been offset by a factor of 
100 for clarity. 

a) the I, Debye scattered intensity, the TE polarized 
ray, and the TE polarized surface wave; 

b) the I 2 Debye scattered amplitude, the TM polarized ray, 
and the TM polarized surface wave. 
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p = 3 intensity near 0 C for the dominant polarization (i.e. the S, 
amplitude) could even be used to experimentally measure the Fock 
transition effect. In Fig. lib, the vanishing of the Fresnel 
reflection coefficient of the geometrical ray at the Brewster angle 
(0“ = 53.12°, 0 B = 115.02°) is also evident. 

A feature of Fig. lib that was unanticipated is the virtual 
nonexistence of supernumeraries for the second order rainbow in the 
nondominant polarization (i.e. the S 2 amplitude) for x - 1000. For 
very large spheres (i.e. x - 20,000) they do appear but are quite 
weak compared to the supernumeraries for the dominant polarization. 
This near-absence of nondominant polarization supernumeraries was 
found for all the rainbows that we examined (the second order through 
fifth order) and seems not to have been commented upon before. The 
absence of the supernumeraries is due to the fact that the Brewster 
angle for the TM polarization is within or near the range of angles 
0, of the supernumerary ray incident closer to the centerline of the 
sphere. The Fresnel reflection coefficient term [R u ] p ' 1 of this ray 
is near zero in much or all of the supernumerary region. Thus the 
ray has a very small amplitude and its interference with the other 
supernumerary ray is quite weak. 

For a plane wave incident on a spherical water droplet in the TE 
polarization state, the intensity of the surface wave - geometrical 
ray interference is about an order magnitude weaker than the 
intensity due to reflection by the droplet surface. For the TM 
polarization state, it is about two orders of magnitude weaker than 
the reflection component. As a result, it would seem that the 
surface wave - geometric ray interference of Fig. 11 should not be 
observable. However, a possibility of observation exists in two 
different experimental situations. First, it has been found that the 
intensity of the second order rainbow for a prolate spheroidal water 
droplet is enhanced with respect to its intensity for a spherical 
droplet. 38 This relative enhancement may render the surface wave - 
geometrical ray intensity comparable to the reflection intensity. 
Second, if a water droplet is illuminated by a laser beam whose 
diameter is smaller than the droplet diameter, shifting the beam off- 
center increases the intensity of the second order rainbow with 
respect to the competing reflected light. 12 For example, consider the 
scattered light at the angle where geometrical optics predicts the 
end of the supernumerary region (i.e. at 0 C = 111.64°). The 
contributing supernumerary rays are incident on the droplet at a 
distance of 0.76a and 1.00a from the center line of the droplet, 
while the rays that are reflected from the droplet are incident 
somewhat closer to the center line of the droplet at a distance of 
0.56a. If the center of an incident laser beam is positioned near 
the edge of a water droplet, the rays that produce the second order 
rainbow will have greater intensity than the rays that are reflected 
at that angle. We experimentally investigate this technique in 
section 8. 

IS 
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Complex Ray - Geometrical Ray Interface 

Consider light rays that make three internal reflections (i.e. 
p = 4 rays) within the sphere as in Fig. 12. The scattering angle of 
the third order rainbow is 0 R = 41.74°, and in ray optics, the 
supernumeraries occur for scattering angles between 0 R and 
0 C = 28.86°. The complex ray contributes to the p = 4 intensity for 
8 > 0 R . In addition, a third geometrical light ray that enters on 
the opposite side of the centerline with 8, < 35° and makes three 
internal reflections, also exits at 0 > 0 R and interferes with the 
complex ray as in Fig. 12. 


The intensity corresponding to the S, scattering amplitude p = 4 
term of the Debye series, shown in Fig. 13, exhibits the third order 
rainbow, a few broad supernumeraries for 0 < 6 R , and a much finer 
oscillatory structure for 0 > 6 R . The amplitude of the finer 
oscillatory structure grows with increasing 6 until the scattering 
angle is about 52°. Thereafter the amplitude decreases. Such a 
behavior is consistent with the interference between the p = 4 
complex ray and the above-mentioned third p = 4 geometrical ray, with 
the complex ray dominating for 0 < 52° and the geometrical ray 
dominating for 6 > 52°. In order to verify this, we compared the 
p = 4 Debye series intensity with 


■i approximate l^third geometrical ray + Airy (6) 


(27) 


in Fig. 13. Again, the comparison of the amplitude and period of the 
oscillatory structure is quite good. 


As to the potential observability of this effect, for scattering 
by a sphere in this angular region the scattering is dominated by 
transmission and reflection. Although the third order rainbow from 
a prolate spheroidal water droplet has been tentatively observed, 3 * 
the complex ray - geometrical ray interference pattern of Fig. 13 is 
more than three orders of magnitude weaker than the rainbow maximum. 
Further, the two rays that contribute to this interference pattern 
enter on opposite sides of the droplet and the dominant p = 0 
reflected ray is incident on the sphere at virtually the same 
location where the third order rainbow ray enters. This renders the 
off-center laser beam technique ineffective in this situation. 
Possibly the only hope for the observation of the complex 
ray - geometrical ray interference would be to use the TM 
polarization state, where the intensity of the reflected ray falls by 
over an order of magnitude while the intensity of the third order 
rainbow falls by about a factor of four. 


Surface Wave - Complex Ray Interface 

Since no geometrical light rays make one internal reflection 
within the sphere and exit with a scattering angle less than 
6 R = 137.92°, the p = 2 term of the Debye series should be a near- 
ideal situation in which to numerically examine weak scattering 
mechanisms in the region 8 < 8 R . Near the first order rainbow, the 
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FIGURE 12: Ap = 4 geometrical ray (solid line) and the p = 4 complex 
ray (dotted line) that interfere in the far-zone. The 
dashed line marked R 3 is the third order rainbow ray. The 
complex ray occurs for scattering angles larger than that 
of the rainbow ray. 
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FIGURE 13: Comparison of the p = 4 Debye component intensity I, with 
the approximation of eq. (27) for the TE polarization state 
and for x = 600 and n = 1.333. The approximation of 
eq. (27) has been offset by a factor of 100 for clarity. 
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p = 2 intensity is dominated by the complex ray which falls off 
faster than exponentially as a function of 6. Since the intensity of 
the secondary radiation shed by surface waves falls off exponentially 
as a function of 6, it has been previously noted that the scattered 
intensity due to surface waves produced on the other side of the 
droplet, as shown in Fig. 14, will eventually overtake the complex 
ray scattered intensity , 45 ^ 6 As a result, an interference pattern will 
be formed in the angular interval in which these two effects are of 
comparable strength. This interference in the p = 2 Debye intensity 
corresponding to the S, scattering amplitude is shown in Fig. 15. For 
0 < 115°, log I, is linear in scattering angle indicating surface wave 
dominance, and for d > 115° log falls off faster than linearly in 
scattering angle indicating complex ray dominance. Again, in order 
to verify this, the intensity of the p = 2 term of the Debye series 
is compared to 

approx ima te ( 6 ) = l^a^a ( 6 ) + E& y ( 0)| 2 (28) 

in Fig. 15. The agreement between eq. (28) and the Debye series 
result is quite good considering that the Airy theory 
parameterization of eq.(9) and the surface wave parameterization of 
eq.(21) are expected to be only qualitatively accurate for size 
parameters as low as x = 100. Since the surface wave-complex ray 
interference structure of Fig. 15 occurs at an intensity level about 
seven orders of magnitude below the intensity of the rainbow peak, 
the effect is not expected to be observable. 

Measurement of Surface Wave - Geometrical Ray Interface 

An experiment was performed to observe the interference between 
geometrical rays and the secondary radiation shed by surface waves. 
The region showing most promise for making this measurement is at 
scattering angles between 90° and 110° in the vicinity of the second 
order rainbow. A focused laser beam was employed as shown in Fig. 16 
and described in Section 5 to increase the visibility of the p = 3 
rainbow and its supernumeraries, and decrease the background 
intensity of the light reflected by the droplet. 

A focused, linearly polarized beam from a five watt argon ion 
laser (\=514.5 jum) was used to illuminate water droplets that were 
generated at a rate of 53,240 droplets per second by a vibrating 
orifice droplet generator as shown in the experimental setup in 
Fig. 17. The size of the droplets was determined by weighing a 
sample of the droplets collected over a period of several minutes. 
Details of this method are given in Ref. 47. Uniform size droplets 
(86.6 nm ± 1.5 jn m) passed through the waist of the laser beam that 
was measured to be 40 jum ± 5 n m . A 25 cm focal length lens was found 
to be ideal for this application for two reasons: 1) it produced a 
waist that was roughly half the diameter of the droplet and 2) the 
length of the waist was long enough to insure that the rays incident 
on the droplet were parallel. If the rays been converging or 
diverging, then the positions of the rainbow and its supernumeraries 
would have been shifted. When a viewing screen was placed beside the 
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FIGURE 14: A p = 2 surface wave (solid line) and the p = 2 complex 
ray (dotted line) that interfere in the far-zone. The 
dashed line marked R 1 is the first order rainbow ray. The 
complex ray occurs for scattering angles smaller than that 
of the rainbow ray. 
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FIGURE 15: Comparison of the p = 2 Debye component intensity I, with 
the approximation of eq. (28) for the TE polarized surface 
wave and complex ray, for x = 100 and n = 1.333. The 
approximation of eq. (28) has been offset by a factor of 
100 for clarity. 
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FIGURE 16: A laser beam with a small diameter, illuminating a water 
droplet close to one edge will cause the p = 3 rainbow and 
its supernumeraries to be more intense than other rainbows 
and confine the reflected rays to the forward scattering 
hemisphere. 
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FIGURE 17: Experimental apparatus. Light scattered by the droplets 
in the angular range 30° < 6 < 150° is observed on the 
viewing screen. 
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droplet stream, and parallel to the laser beam, the scattering 
interval 30° < 9 < 150° could be observed, and the first and second 
order rainbows were readily evident. By adjusting the droplet stream 
with the use of a micropositioner so that the laser beam was incident 
on either the left or right side of the droplet centerline, either 
the first or second order rainbow could be made more intense while 
the other rainbow was extinguished. When the laser beam was incident 
near the side of the droplet closest to the screen, thus illuminating 
the second order rainbow, the reflected intensity was largely 
confined to the forward hemisphere and no longer obscured the second 
rainbow supernumeraries. A photograph of the second order rainbow 
using this arrangement is shown in Fig. 18. The photograph was taken 
by replacing the screen with a film holder. 

A frame of polaroid film was exposed by the scattered light for 
0.005 sec by means of a camera shutter located at the exit port of 
the laser. Longer exposures using lower power laser beams were not 
taken in order to avoid blurring produced by small fluctuations in 
the droplet stream from air currents in the room. Moving the droplet 
generator head very close to the laser beam to reduce the effects of 
the air currents was not practical because the droplets were unstable 
and not spherical in this region. For the conditions of Fig. 18 the 
droplet generator head was approximately 5 cm from the laser beam. 
Note that the figure is oriented so that the scattering pattern is 
consistent with the orientation of the rays in Figs. 16 and 17, thus 
the scattering angle in Fig. 18 increases from right to left. In 
Fig. 18 the broad intense region is the p = 3 rainbow, with the 
supernumerary region located to the right of it in the direction of 
decreasing scattering angle. Farther into the forward scattering 
direction, the supernumeraries become washed out by the increasing 
background of the p = 0 reflected light. 

After photographing the scattering pattern, the photograph was 
digitized and compared with the Debye series calculation for p = 3 as 
shown in Fig. 19. No compensation for the nonlinearity of the film 
was made. The abscissa in this figure has two axis labels: one in 
millimeters and one in degrees. A scattering angle of 90° 
corresponds to a position on the film of 0 mm. The degree tick marks 
are not linear because the film was held with a flat rather than a 
circular-shaped film holder. The experimental data shows 

interference oscillations which begin at the rainbow angle, continue 
through the two-ray supernumerary region and extend well into the 
one-ray region which begins on the film at the position -12 mm. The 
positions of these oscillations correspond to the expected positions 
of the supernumeraries of the second order rainbow for an x = 528.8 
sphere as obtained from the p = 3 Debye component intensity. 
Additional calculations were performed to be certain these 
oscillations were not a result of interference between the one 
remaining p = 3 geometrical ray and the p = 0 reflected rays. These 
calculations showed that such an interference pattern would have had 
a greatly different periodicity than the one observed. Thus we are 
confident that the oscillations that were observed in the one-ray 
region are a result of the interference between the geometrical ray 
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FIGURE 18: Scattered light intensity for a 40 /xm diameter laser beam 
incident near the edge of an 86.6 /xm diameter water 
droplet. The scattering angle increases from right to 
left. The broad illumination on the right side of the 
photo is reflection from the droplet. The interference 
pattern on the left side of the photo is the second order 
rainbow and its supernumeraries. 
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FIGURE 19: Comparison between the Debye component intensity and 
experiment. 

a) Detector response (normalized) as a function of 
position on the film for the digitation of the 
photograph of Fig. 18. The four oscillations on 
the film between -12 mm and -6 mm are the geometrical 
ray - surface wave interference. The rise in the 
detector response between -6 mm and 0 mm is due to 
reflection by the droplet surface. 

b) The p = 3 Debye intensity I, for x = 528.8 and 
n = 1.333 as a function of scattering angle. 
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and the secondary radiation shed by the surface waves described in 
section 5. 


Conclusions 

When one thinks about electromagnetic scattering by a dielectric 
sphere, one's physical picture of the scattering is very different 
than one's mathematical calculation of it. Physically, one thinks of 
scattering as the resultant sum of a number of "mechanisms" such as 
geometrical rays, rainbows, glories, surface waves, and diffraction, 
and the wave interference between them. Mathematically, one 
calculates the scattering in terms of the effect of the interface 
between two dielectrics on spherical multipole partial waves and the 
interference between the partial waves. The difficulty in 
reconciling these two modes of thought lies in the fact that 1) many 
partial waves contribute to a single physical mechanism and 2) only 
a portion of the scattering amplitude for each partial wave 
contributes to a given mechanism. The connection between the 
mathematical and physical modes of thought is expressed in a simple 
way by the Debye series. We hope to have demonstrated the simplicity 
and power of this connection in this paper. 

The Debye series also makes evident one other feature of 
electromagnetic scattering which is not evident in the geometrical 
model of scattering, namely the scattered intensity as a function of 
scattering angle is smoothly varying. Ray theory allows (in fact 
demands) discontinuities in both the scattered intensity as a 
function of scattering angle and its derivative. The discontinuity 
in the intensity occurs at the two-ray to zero-ray transition at the 
rainbow scattering angle. The discontinuity in the derivative of the 
intensity occurs at the one-ray to zero-ray transition at the 
critical scattering angle of the grazing incident ray. Wave theory 
smooths these discontinuities by the action of the complex ray at the 
two-ray to zero-ray transition and the secondary radiation shed by 
surface waves at the one-ray to zero-ray transition. These physical 
mechanisms which implement the smoothing of the scattered intensity 
at the more-rays to fewer-rays transitions produce new features in 
the scattering as well. They interfere with each other and with 
geometrical rays, giving additional structure to the scattered 
intensity. Under favorable circumstances, at least one of these new 
structures, the geometrical ray - surface wave interference is 
observable in laser light scattering. 
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Appendix A - Derivation of the Debye Series 

Consider monochromatic transverse electric (TE) and transverse 
magnetic (TM) spherical multipole waves in a medium whose refractive 
index is n. The time dependence of the waves is taken to be e''“ l . The 
electric and magnetic fields of the spherical waves are derived from 
a scalar potential \p(r,0,(p) which satisfies the wave eguation 

V 2 v|/ + n 2 k 2 \\i = 0 (Al) 

where 

k = — . ( A2 ) 

c 

The fields are 
E = -r x v\|j 

B = — — v x E, 

0) 

for the TE spherical multipole waves, and are 

— / c? — 

E = — — v x B 
(0 n 2 

B = — r x v\ Jr, 

c 


(A3) 


(A4) 


for TM spherical multipole waves. The solutions to the scalar wave 
eguation (Al) are 


i|r (r, 0, <J>) 


= E 


J J*i (nkr) 1 


P”(cos0) 


COS 7T2 4> 

sin /n 4> 


(A5) 


where j f and n f are spherical Bessel functions and spherical Neumann 
functions respectively (also known as spherical Bessel functions of 
the first and second kind) and P, m are associated Legendre 
polynomials. We adopt the sign convention for the spherical Bessel 
and Neumann functions 


j Q ( nkr ) 


sin ( nkr ) 

nkr 


n 0 (nkr) 


-cos (nkr) 

nkr 


(A6) 


With this sign convention, the Wronskian relation for the spherical 
Bessel and Neumann functions is 


j,(nkr) n[(nkr) - j' t (nkr) n,(nkr) = (nkr) 2 , ( A7 ) 

where the symbol prime indicates the derivative of a function with 
respect to its argument. The solutions of eq. (A5) represent standing 
waves. When considering outgoing and incoming traveling waves, we 
employ spherical Hankel functions of the first and second kind 


jr 
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h, (nkr) = j t (nkr) + i n t (nkr) 

(A8) 

h, u) (nkr) = j t (nkr) - i n t (nkr) 
respectively. 

Consider a dielectric sphere whose radius is a and refractive 
index is n, (region 1) embedded in another dielectric material whose 
refractive index is n 2 (region 2) . The center of the sphere is taken 
to be at the origin of the coordinates. If a spherical multipole 
wave propagates from one region to another, the tangential components 


of its electric and magnetic 
interface. For TE waves, this 
requires that 

Y(x) = Y(y) 

n 2 

and 

Y'U) = Y 7 (y) 
where 
x = n 2 ka 

y = n^ka 
and 

Y (nkr) = nkr \Jj (nkr) . 

For TM waves, the continuity of 
YU) = Y (y) 
and 

Y'(x) = Y ; (y) 
n 2 n x 


fields must be continuous at the 
continuity of the field components 

( A9 ) 

(A10) 

(All) 

(A12) 

the field components requires that 

( A13 ) 

(A14) 


Consider the single incoming TE or TM spherical multipole wave 
in region 2 

(cos /r?4>) 

Y = H t {2) (n 2 kr) P, m ( cos0W (A15) 

[sin/n4>J 

where 

tf/ 2) (n 2 kr) = n 2 kr h, U) (nkr) . (A16) 

When this spherically incoming wave encounters the interface at 
r = a, a portion of it, T f 21 , is transmitted into the sphere and a 
portion of it, R, 22 is reflected back into region 2. The complete TE 
or TM' multipole wave in the two regions is then 
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T, = if 1 H, U) (n x kr) P, m (cos0) | cos/7?(t> l 

[sin/n4> J 

T 2 = [h, {2) ( n 2 kr ) + Rf 2 H t ll) ( n 2 kr )] 


for r <. a 


(A17) 


/or r £ a . 


„ (cos/n<t>l 

X P, (COS0) { i 

[sin7i?4> J 

Applying the boundary conditions (A9-A10) or (A13-A14) yields two 
linear equations in the two unknowns T, 21 and R, 22 . The solutions of 
these equations are 


( A18 ) 


T, f l = 

-fM 

2i 

\ n 2) 

D 

and 



R 22 = 

[a/// 2)/ 

(x) H t (2) (y) - P H t {2) (x) 


where 


a = 


1 

no 


D 

for TE spherical multipole waves 
for TM spherical multipole waves f' 


(A19) 


(A20) 


P - 


and 


£ i 

n o 


for TE spherical multipole waves 
for TM spherical multipole waves 


(A21) 


(A22 ) 


D = -a H t n) '(x) H t {2) (y) + p H, ll> (x) H t {2) '(y) , 

\ 

and where eq. (A7) has been used to simplify the expression for T, 21 . 

Consider now the single outgoing TE or TM spherical multipole 
wave in region 1 


Y = H, U) (n r kr) P f m (cos6) 


[cos m<J>] 
[sin/r)4>j 


(A23) 


When this spherically outgoing wave encounters the interface at 
r = a, a portion of it, T, 12 , is transmitted out of the sphere and a 
portion of it, R f n , is reflected back into it. The complete TE or TM 
multipole wave in the two regions is then 
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Y, = [O (n.kr) + R, 11 H t i2) (n.kr)] P, m (cos0) | cos/n< ^). for 
1 J (sin /77 4>J 


r <. a 


(A24) 


T 2 = Tf 2 H, U) (n 2 kr) P t m (c os0) 


( 003/774)1 

[sin/774) 


f 


for r £ a. 


Applying the boundary conditions (A9-A10) or (A13-A14) again 
yields two linear equations in the two unknowns T, 12 and R," . The 
solutions of these equations are 


if 2 = 
and 

-2 i/D 

(A25) 

l?. 11 = 

[ ocH^'ix ) H, ll) (y) - p h} 1] (x) H, U) '(y)] 

(A26) 


D 

Eqs . ( A18-A22 , A25-A26) are identical to eqs.(2.11, 2.13, 2.15-2.18) 
of Ref. 8 except that the right hand side of eqs . (2 . 17-2 . 18 ) of that 
reference should be multiplied by nx 2 . 

Consider the quantity ( 1-R, 11 ) ( 1-R f 22 ) -T f 21 T f 12 for either the TE or 
TM spherical multipole waves. Substituting eqs.(A19, A26) , and using 
the expressions for T, 21 and T/ 2 before the simplifications provided by 
the Wronskian relation, the quantity is found after much algebra to 
factor into 


(1-0 (1-0 ~ TfO = 4 f + P«7 t (x) j/(y)| 

where 

J t (nkr) = nkr j ( (nkr) . 

Dividing eq. (A27) by (1-R f n ), we obtain 


( A27 ) 


(A28 ) 


-oc jj (x) J t (y) +p J, (x) jj (y) 

-a Hf 1) ' {x) J ( [y) +P H t {1) ( x ) Jiiy) 


1 - R , 22 - 


Tf 1 !? 2 


1 - 1 ?, 


11 


(A29) 


1 

2 


1 -R 22 - £ if 1 (0*- 1 Tt‘ 

p- 1 


The left hand side of eq. (A29) is recognized as the Mie far- 
field partial wave scattering amplitudes a ( (TM amplitude) and 
b f (TE amplitude) . The Mie scattered electric field is written in 

terms of the a, and b, partial wave amplitudes as 

» 
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(A30) 


E <3 

Z?(i?,0,(}>) = e iw [-i5 2 (0) cos<t> u Q + i S 1 (0) sin<{) u^] 


for an incident electric field in the G* direction where 




ft 

(A3 1 ) 

s 2 ( 0 ) = «■ We>>, 


( A32 ) 

n,(0) = 1 -Pj 1 ( 0 ) , 

* sin0 


(A33) 

t,(6) " is p,1(e) ' 


(A34) 

and 



lax = x + 4.05X 1 / 3 + 2, 


(A35) 

also 



1,(0) = |5, (0) | 2 , 


(A36) 

J 2 (0) = |S 2 (0) | 2 . 


(A37 ) 

For a large sphere x >> 1, large partial waves Z » 1, and 6 away 
from 0° and 180°, the angular function n t (6) is small in comparison 
to T e ( 0). 48 In this limit, the scattering amplitude S, becomes 
associated with the TE spherical multipole waves and the TE polarized 
geometrical rays, while the scattering amplitude S 2 becomes associated 
with the TM spherical multipole waves and the TM polarized 
geometrical rays. 


The right hand side of eq. (A29) is the Debye series expansion of 
the partial wave scattering amplitudes. The various terms of the 
right hand side of eq. (A29) have simple physical interpretations. 
The first term, %[1], when inserted for a f and b, in the expressions 
for the Mie electric field and summed over Z, describes the 
diffraction of the incident plane wave around the sphere. 49 

The second term, %[-R f 22 ], when inserted for a t and b, and summed 
over Z, represents the outgoing spherical multipole waves that have 
reflected from the surface of the sphere. The diffraction and 
reflection terms when taken separately, do not approach zero as Z 
goes to infinity (and thus the sums over Z go to infinity) . However, 
if the two terms are added together then summed over Z they do 
approach zero, which leads to a finite sum over Z. 3 ‘ Thus, in its 
present form the Debye series is not amenable to separation of the 
diffraction term from the reflection term. Previous calculations 
have been made in which the diffraction term was summed separately 
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from the reflection term 49 ' 50 by truncating the sum over Z at Z = x. 
Truncation in this manner is an approximation and is only accurate 
for very large values of x. Truncation of the diffraction term of 
Debye series at Z = x has been shown to be equivalent to Fraunhofer 
diffraction by a circular aperture of radius a. We conjecture that 
the separate divergences of the diffraction and the p = 0 terms of 
the Debye series are artificial in that they result from the infinite 
extent of the incident plane wave. We believe that for an incident 
beam of finite lateral extent 51 (e.g. a Gaussian) the sums over Z 
should converge individually. 

The third term of eq. (A29) is an infinite sum which describes a 
geometric series. Individual terms of the series represent light 
that has penetrated the sphere, undergone p-1 internal reflections, 
then exited the sphere. For example, when calculating the p = 2 
contribution to the Debye series, the appropriate term in the sum 
that describes a, and b, is %[T, 21 R, 11 T, 12 ] . This represents the 
spherical waves that were transmitted into the sphere [T, 21 ], 
propagated radially inward toward the center of the sphere, passed 
through it, propagated radially outward, reflected off the inner 
surface of the sphere [R f “], propagated inward toward the center and 
back out a second time, and finally transmitted out of the 
sphere [T f 12 ] and propagated into the far-zone. 

If all the Debye series terms are first added together then 
summed over Z, the results are identical to those from Mie 
scattering. In this sense, eq. (A29) interprets the Mie partial wave 
scattering amplitudes as the composite result of a multiple 
scattering process. 
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